Observation of double radiative capture on pionic hydrogen. 
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We report the first observation of double radiative cap- 
ture on pionic hydrogen. The experiment was conducted at 
the TRIUMF cyclotron using the RMC spectrometer, and de- 
tected 7-ray coincidences following 7r~ stops in liquid hydro- 
gen. We found the branching ratio for double radiative cap- 
ture to be (3.05 ± 0.27(stat.) ± 0.3I(syst.)) x 10~ 5 . The mea- 
sured branching ratio and angle-energy distributions support 
the theoretical prediction of a dominant contribution from the 
7T7T — > 77 annihilation mechanism. 
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Negative pions stopped in hydrogen form pionic hy- 
drogen atoms. These atoms can disintegrate via several 
modes that include the well-known processes of charge 
exchange 7T~p — > 7T°n [1], radiative capture 7r~p — > 771 
[1], and pair production tt~p — > e + e~n [2,3]. 

However, for pionic hydrogen an additional mode of 
capture is predicted by theory, 



7T p — * 777I 



(1) 



This double-radiative process has been investigated theo- 
retically by several authors, including Ericson and Wilkin 
[4], Christillin and Ericson [5], Gil and Oset [6], and 
Beder [7]. The predicted branching ratio is 5.1 x 1CT 5 
[7], with a mechanism that is dominated by the annihi- 
lation of the stopped, real n~ on a soft, virtual 7T + , i.e. 
7T7T — > 77. Beder also predicted different photon energy- 
angle distributions for the contributing annihilation and 
bremsstrahlung mechanisms. 

The underlying dynamics of 7T7t annihilation in double 
radiative capture is rather intriguing. For example, it led 
Ericson and Wilkin [4] to suggest the reaction as a probe 
of the pion field in the nucleus, and Gil and Oset [6] to 
suggest the reaction as a novel window on the 7T7t — > 77 
vertex. Also, the related "/p —* 77m reaction was consid- 
ered by Wolfe et al. [9] and Drechsel and Fil'kov [10] as 
a possible probe of the pion polarizability. 

The only experimental search for double radiative cap- 
ture on pionic hydrogen was conducted by Vasilevsky et 
al. [11] at JINR. They used a large-acceptance photon- 
pair spectrometer and obtained a branching ratio upper 
limit of 5.5 x 10~ 4 . However, double radiative capture on 
beryllium and carbon has been observed in experiments 
by Deutsch et al. [12] at CERN and Mazzucato et al. [13] 



at TRIUMF. Unfortunately, these data are difficult to in- 
terpret due to i) nuclear structure effects and ii) capture 
occurring from both the s- and p-states of the 7rBe and 
7rC atoms. 

Our experiment was performed at the TRIUMF cy- 
clotron using the RMC spectrometer [14]. The incident 
beam had a pion flux of 7 x 10 5 s _1 , a central momen- 
tum of 81.5 MeV/c, and electron and muon contamina- 
tion of 18% and 9% respectively. The incoming pions 
were counted in a 4-element plastic scintillator telescope 
and stopped in a 2.7 liter liquid hydrogen target of length 
15 cm, diameter 16 cm, and wall thickness 254 /im [14]. 
The outgoing photons were detected by pair production 
in a 1 mm-thick cylindrical Pb converter and electron- 
positron tracking in cylindrical multiwire and drift cham- 
bers. A 1.2 kG axial magnetic field was used for momen- 
tum analysis and concentric rings of segmented scintilla- 
tors were used for fast-triggering. The trigger scintilla- 
tors comprised the A-ring (just inside the Pb converter 
radius), the C-ring (just inside the multiwire chamber 
radius), and the D-ring (just outside the drift chamber 
radius). For more information on the RMC spectrometer 
see Wright et al. [14]. Note that in this experiment we 
moved the Pb converter from just inside the C-counter 
radius to just outside the A-counter radius. 

For 7r~p — > 7771 data-taking we employed a two- 
photon trigger based on the hit multiplicities and the hit 
topologies in the trigger scintillator rings and the drift 
chamber cells. A typical 7T~p — > 7771 event that fulfilled 
the trigger is shown in Fig. 1. It has zero hits in the 
A-counter ring, two hits in the C-counter ring, and four 
hits in the D-counter ring. To reduce the high rate of 
back-to-back photons from ir° — > 77 decay, we rejected 
photon-pairs reconstructed with drift cell hits or trigger 
scintillator hits separated by large azimuthal angles. 

During a four week running period we collected ir~p — > 
77n data from a total of 3.1 x 10 11 pion stops in liquid 
hydrogen. Calibration data with a dedicated n° — » 77 
trigger were also taken periodically. 

One source of background was real 7-7 coincidences 
arising from tt° — > 77 decay. The 7T°'s were produced 
by either at-rest or in-flight pion charge exchange. The 
at- rest source yields 7T°'s with energy T = 2.8 MeV and 
decay photons with opening angles cos (9 < —0.91, while 
the in-flight source yields 7r°'s with T < 15 MeV and 
photons with cos 9 < —0.76. The at-rest background was 
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roughly 1600 x the double-radiative capture signal, and 
the in-flight background was roughly 10 x the double ra- 
diative capture signal. Consequently, for n~p — ► -fjn the 
cos 9 < —0.76 region was overwhelmed by tt° background, 
and due to the finite resolution of the photon-pair spec- 
trometer, the 7r° background was a potential problem for 
opening angles with cos9 > —0.76. 




FIG. 1. A typical n~p — * 7771 event. The plot shows the 
fit in the plane perpendicular to the beam axis. The elec- 
tron-positron pairs converge at the lead converter and the 
reconstructed photon pairs originate from the hydrogen tar- 
get located at the center. The trigger pattern of zero hits in 
the A-counter ring, two hits in the C-counter ring, and four 
hits in the D-counter ring is also displayed. For scale, the 
radius of the D ring is about 60 cm. 

Another source of background was accidental 7-7 co- 
incidences arising from simultaneous multiple 7r~ stops. 
The pion beam had a micro-structure with a pulse width 
of 2-4 ns and a pulse separation of 43 ns. With an in- 
cident flux of 7 x 10 5 s" 1 the probability for more than 
one pion arriving in a single beam pulse is 1.5%. Multi- 
ple pion stops in one beam pulse can yield a 7-ray pair 
by the accidental coincidence of one photon from each 
pion. This background was roughly 900 x our signal. 
It yields photon-pairs with opening angles 0-180° and 
summed energies 106-258 MeV. Note that the summed 
energy from random background events can exceed the 
m^-limit for single tt capture. 

In analyzing the data a number of cuts were applied 
to identify photon pairs and reject background sources. 
A tracking cut imposed minimum values for the number 
of points in the tracks and maximum values for the chi- 
squared of fits to the tracks. A photon cut required that 
the electron-positron pairs intersect at the Pb converter 
and that the reconstructed photon pairs originate from 
the H2 target. To reject the multi-7r background we im- 
posed a C-counter timing cut and a beam telescope am- 
plitude cut. The telescope amplitude cut was imposed 
on the normalized sum of the light output from the eight 
photo-multiplier tubes viewing the four individual beam 
scintillators. The ±4 ns C-counter timing cut was im- 
posed on the time difference between the two C-counters 
intersecting the two emerging e + e~ pairs. The remain- 
ing inefficiency in rejecting the accidental coincidences 



was 1.3 x 10~ 4 . To reject the tt° background we imposed 
a photon opening angle cut of cos 9 > —0.1. 

A total of 2.3 x 10 6 photon pairs passed both the 
tracking cuts and photon cuts. These photon pairs are 
shown in Fig. 2 and are dominated by the backgrounds 
from 7r° decays and multi-7r stops. The multi-7r back- 
ground is clearly seen in the summed energy spectrum 
as events with E sum > 150 MeV and the n° background 
is clearly seen in the opening angle spectrum as events 
with cos 9 < —0.76. The beam telescope amplitude cut 
removed about 0.8 x 10 6 accidental 7-7 coincidences from 
multi-7T stops, and the photon opening angle cut removed 
about 1.4 x 10 6 real 7-7 coincidences from ir° decays. A 
total of 635 events with E sum > 80 MeV and cos 9 > —0.1 
were found to survive all cuts (see Fig. 3). 
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FIG. 2. The summed energy spectrum (top) and open- 
ing angle spectrum (bottom) for reconstructed photon-pairs, 
i.e. the events passing the tracking and photon cuts. 
Note that the multi-7r background can produce events with 
E 3um > 150 MeV (see top plot) and the n° background will 
produce events with cos 6 < —0.76 (see bottom plot). The 
Monte Carlo generated 7r° background is shown overlayed as 
the shaded histogram in the bottom plot. The arrow in the 
upper plot indicates the E sum > 80 MeV cut and the arrow 
in the lower plot indicates the cos 6 > — 0.1 cut. 

A small quantity of two-photon background from n° 
decays and from multi-7r stops does however survive the 
applied cuts. The remaining ir° contamination was sub- 
tracted using (i) the observed number of n° events with 
cos 9 < —0.76 and (ii) the known angular response of the 
photon-pair spectrometer. The remaining multi-7r con- 
tamination was subtracted using (i) the observed num- 
ber of 27T events with E sum > 170 MeV and (ii) the 
measured sum energy spectrum for the multi-7r back- 
ground. These procedures indicated 53±30 7r° back- 
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ground events, or (8.3±4.8)%, and 100±16 multi-7r back- 
ground events, or (15.7 ± 2.5)%, with E sum > 80 MeV 
and cos 9 > —0.1. After subtraction this yielded a to- 
tal of 482±42 7r~p — ► 7777 events with E sum > 80 MeV 
and cos# > —0.1. Using the data from Refs. [12,13] we 
further estimated the backgrounds originating from the 
nuclear (it, 2j) reaction on the target walls, etc, to be 
< 1%. 

In order to compare the experiment with theory we 
performed both measurements and simulations of the 
two-photon response function of the RMC spectrome- 
ter. To measure the response we employed the multi- 
7r accidentals. Specifically, by counting the numbers of 
incoming multi-7r stops and outgoing 7-7 accidentals, we 
mapped the detector's response versus energy and angle 1 . 
Note that the calibration data and 7r~p — > 7771 data were 
collected simultaneously and passed through the same 
cuts (with the exception of the beam telescope cut). Con- 
veniently, the energy range and angular range for multi-7r 
accidentals covers the kinematical range for ir~p — ► 7771 
events. 

To simulate the detector response function we used a 
Monte Carlo program. The program incorporated both 
the detailed geometry of the RMC detector and the de- 
tailed interactions of the relevant particles. Our program 
was based on the CERN GEANT3 package [16] and is de- 
scribed in Wright et al. [14]. We tested the simulation by 
comparison to the response function measurements with 
the multi— 7r accidentals. We found the energy-angle dis- 
tributions from experiment and simulation to be in very 
good agreement. However the absolute detection efficien- 
cies from experiment and simulation were found to differ 
on average by 10%. Moreover, the measured acceptance 
was found to vary by ±4% from run to run and decreased 
by 10% between the early runs and the late runs. We 
attributed these variations to changes in the chamber ef- 
ficiencies and fluctuations in the chamber noise, neither 
of which were incorporated in the simulation of the ac- 
ceptance. To account for these differences between the 
measured and the simulated acceptance, we employed a 
multiplicative correction factor of F = 0.90 ± 0.09. The 
quoted uncertainty is very conservative, and embodies 
the entire variation of the measured acceptance over the 
running period. 

In addition we compared the results of measurements 
and simulations of photon-pairs from at-rest charge ex- 
change ir~p — > 7r°77 followed by neutral pion decay 



1 A minor complication is that the two photons from the 
multi-7r stops may have time differences of up to 4 ns (i.e 
the pion beam pulse width). This time difference slightly 
decreases both the track reconstruction efficiency and the two- 
photon acceptance. From simulations we found the loss in 
acceptance to be 6%. 



7r° — ► 77. Here we used dedicated l ir° runs' with a mod- 
ified trigger arrangement for these back-to-back photon 
pairs. We found excellent agreement between the simu- 
lation and the measurements using the factor F = 0.90. 
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FIG. 3. Comparison of the opening angle distributions 
from the background subtracted experimental data (open cir- 
cles) and the theoretical calculation (curves). The dashed 
curve is the -kit annihilation process, the dotted curve is the 
NN bremsstrahlung process, and the solid curve is the full 
calculation. These curves are convoluted with the response 
function of the RMC spectrometer. 

The branching ratio for double radiative capture on 
pionic hydrogen was obtained via 



N n - ■ eft ■ F ■ Cb m ■ c s top 

where N n - is the number of livetime-corrected pion 
stops, is the number of background-subtracted 

tt~P — > 7777 events, and efl ■ F is the detector accep- 
tance. Note that the appropriate acceptance was ob- 
tained using Monte Carlo [15] with the 7r~p — > 7777 kine- 
matical distributions taken from Bcdcr [7]. The factor 
Cstop = 0.85 ± 0.01 accounts for the fraction of incident 
pions that stopped in hydrogen (see Wright et al. [15] 
for details) and the factor C(, TO = 0.99 accounts for the 
efficiency of ir~p — > 7777 events passing the beam tele- 
scope cut. Using Eqn. 2 we obtained a branching ratio 
of (3.05±0.27(stat.)±0.31(syst.)) x 10~ 5 . Note that the 
quoted uncertainty contains a statistical error of ±8% 
from A^ 77 and a systematic error of ±10% in total. The 
systematic error is completely dominated by the ±10% 
uncertainty in the determination of the acceptance efl-F. 
The uncertainties in N v - , c st0 p and C(, TO were each < 2 
% and entirely negligible. We stress that the result we 
quote is the total ir~p — ► 7777 branching ratio for all 
photon energies (0 < E 1 < m v ) and all opening angles 

(-1.0 < cos e < +1.0). 

Our quoted branching ratio was extracted assuming 
the energy-angle distributions calculated by Beder [7] , al- 
though we actually only observed the region with cos 9 > 
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—0.1 and E-y > 25 MeV. Tests of the sensitivity of the ex- 
tracted branching ratio to the energy-angle cuts revealed 
only a ±2.9% variation for —0.2 < cos 6 < 0.0 and a 
±2.8% variation for analyses with different sum energy 
cuts. Also, when using a phase space energy-angle distri- 
bution rather than Beder's energy-angle distribution [7], 
the extracted branching ratio changed by only —7%. 

In Beder's calculation [7] of double radiative capture 
the main contributions originate from tttt annihilation 
graphs, NN bremsstrahlung graphs, and their interfer- 
ence. The 7T7T annihilation graphs alone account for 
64% of the total branching ratio and yield a distribu- 
tion that is peaked at small opening angles. The NN 
bremsstrahlung graphs alone account for 20% of the 
total branching ratio and yield a distribution that is 
peaked at large opening angles. At threshold, the pion 
bremsstrahlung contributions vanish, and effects of vec- 
tor meson exchange and delta resonance excitation arc 
calculated to be very small. 

In Fig. 3 we compare our measured data with Beder's 
calculation [7]. The background contributions of (8.3 ± 
4.8) % from tt° decay events and (15.7±2.5)% from multi- 
7r stop events have been subtracted from the measured 
data and the resulting 482 events are plotted as open cir- 
cles. The theoretical curves have been convoluted with 
the response function of the RMC spectrometer. The plot 
shows that the ir~p — > -f-fn branching ratio and opening 
angle distributions from experiment and theory are in 
reasonable agreement. The overall consistency of experi- 
ment and theory supports the theoretical prediction of a 
dominant tttt annihilation mechanism. As seen in Fig. 3, 
the bremsstrahlung graphs alone underpredict the data 
by about a factor of five. 

However, our measured branching ratio is somewhat 
smaller than the theoretical branching ratio. We note 
that Beder's calculation was performed at tree-level and 
neglects contributions from pion loops, etc. We there- 
fore speculate that higher-order terms might explain the 
difference between the experimental and the predicted 
branching ratio. A new calculation of double radiative 
capture on pionic hydrogen using chiral perturbation the- 
ory is currently underway [17]. 

We remind the reader of the results from the previ- 
ous nuclear (77,27) measurements. For 12 C, Deutsch et 
al. [12] obtained a partial branching ratio of (1.4 ±0.2) x 
10~ 5 , for E 1 > 25 MeV and cos6» < 0.71, and Maz- 
zucato et al. [13] obtained a partial branching ratio of 
(1.2 ± 0.2) x 10~ 5 , for E 7 > 17 MeV and cos (9 < 0.71. 
Unfortunately the comparison of the earlier nuclear data 
with our hydrogen data is difficult as i) tt capture is pre- 
dominantly from the IS state in 1 H and the 2P state in 
12 C and ii) the nuclear data were mainly taken at large 
two-photon opening angles and our hydrogen data were 
mainly taken at small opening angles. 

In summary, we have made the first measurement of 
double radiative capture on pionic hydrogen by record- 



ing 7-ray coincidences from tt~ stops in liquid H2. We 
found the branching ratio to be (3.05 ± 0.27(stat.) ± 
0.31 (syst.)) x 10~ 5 by assuming the kinematical distribu- 
tions from Beder [7]. Moreover, the measured branching 
ratio and opening angle distribution support the theo- 
retical hypothesis of a tttt annihilation mechanism. We 
hope this work will stimulate further studies into using 
double radiative capture as a novel probe of the proton's 
pion cloud and the tttt — > 77 vertex. 
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